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Abstract 
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■ We study the scaled variances of multiplicity fluctuations in nucleus-nucleus colli- 
', sions at SPS and RHIC energies within the HSD transport model. The HSD results 

are compared with proton-proton data and with predictions of the hadron-resonance 
gas statistical model. We find that the HSD scaled variances u){ for negatively, pos- 
itively, and all charged hadrons in central nucleus-nucleus collisions remain close to 
the LOi values in proton-proton collisions and increase with collision energy as the cor- 
responding multiplicities per participating nucleon. The statistical model predicts 
very different behavior of ui{. However, a comparison with preliminary NA49 data 
for the most central Pb+Pb collisions at SPS energies does not permit to distinguish 
the HSD and statistical model results. New measurements of the multiplicity fluc- 
tuations in nucleus-nucleus collisions in a wide energy region with large acceptance 
are needed to allow for a proper determination of the underlying dynamics. 
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The event-by-event fluctuations in high energy nucleus- nucleus (A+A) col- 
lisions (see e.g., the reviews [T|l2] ) are expected to be closely related to the 
transition between different phases of the QCD matter. Measuring the fluc- 
tuations one might observe anomalies of the onset of deconfinement J3] and 
dynamical instabilities when the expanding system goes through the 1-st or- 
der transition line between the quark-gluon plasma and the hadron gas [I]. 
Furthermore, the QCD critical point may be signaled by a characteristic pat- 
tern in fluctuations [5j. In the present paper we calculate the multiplicity 
fluctuations in central A+A collisions at SPS and RHIC energies within the 
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Fig. 1. The multiplicity (left) and scaled variance (right) of all charged hadrons in 
p+p inelastic collisions as functions of collision energy. The dashed lines correspond 
to the parametrization ([2]) from Ref. [lj. The solid lines are the HSD results. 

microscopic Hadron- String-Dynamics (HSD) transport model [6J which gives 
rather reliable description for the inclusive spectra of charged hadrons in A+A 
collisions from SIS to RHIC energies [7]. 

For a quantitative measure of the particle number fluctuations it is convenient 
to use the scaled variances, 



L)i = 



(Nj) ~ W 
(Ni) 



where (• • • ) denotes event-by-event averaging and the index % means "+", 
and "ch", i.e negative, positive, and all charged final state hadrons. The com- 
pilation of proton-proton (p+p) data for (N c h) and u c h are taken from Ref. [1] 
and presented in Fig. [U The energy dependence can be parameterized by the 
functions [I]: 
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where ^snn is the center of mass energy. 



The HSD model description of the p+p data (for p+p reaction this is al- 
most equivalent to the Lund-String model [8]) is shown in Fig. [l]by the solid 
lines. It gives a good reproduction of the p+p data for (N c h), but slightly un- 
derestimates uj c h at high collision energies. For negative and positive charged 
hadrons the average multiplicities and scaled variances in p+p collisions can 
be presented in terms of the corresponding quantities for all charged particles, 
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In general, one can define two groups of hadron observables. The first group 
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Fig. 2. The HSD results for n c h {left) and u c h {right) in p+p (dashed line), p+n 
(dashed-dotted line), n+n (dotted line), and N+N (solid line) collisions. 

includes observables which are rather similar in A+A and p+p collisions, 
thus, they can be reasonably described within the wounded nucleon model 
(WNM) [TT] which treats the final state in A+A collisions as the result of 
independent nucleon-nucieon (N+N) collisions. The second group consists of 
A+A observables which are very different from those in p+p collisions (e.g. 
collective flow, multi-strange baryon production). The question arises: Are the 
multiplicity fluctuations in A+A collisions close to those in p+p reactions, or 
are they very different? 

To answer this question let us first consider the model predictions. To compare 
central collisions of heavy nuclei and N+N collisions within the HSD model 
we construct the multiplicities and scaled variances of N+N reactions using 
the HSD results for p+p, p+n and n+n collisions: 

(N» N ) = a pp (NT) + a pn (Nf n ) + a nn (iV™) , (4) 

= [«RP «T ( N * P ) + <*P» «T ( N D + OCnn (NT)] , (5) 

where a pp = Z 2 / A 2 = 0.16, a pn = 2Z{A - Z)/A 2 = 0.48, a nn = {A - 
Z) 2 1 A 2 = 0.36 are the probabilities of proton-proton, proton-neutron, and 
neutron-neutron collisions in Pb+Pb (A=208, Z=82) or Au+Au (A=197, 
Z=79) reactions. The results for p + p,p + n and n + n collisions - presented 
in Fig. [2] - are very close to each other. A small difference between p+p and 
N+N collisions (< 5%) is only seen at SPS energies (shown separately in the 
upper left corners) and gradually disappears at RHIC energies. 

In Fig. [3] the HSD model results are shown for the multiplicities per partici- 
pating nucleons, rij = (Nj) /(N P ), and for the scaled variances, u i: in central 
collisions (zero impact parameter, b = 0) of Pb+Pb at E iab = 10, 20, 30, 40, 
80, 158 AGeV and Au+Au at v /s^7 = 62, 130, 200 GeV. From Fig. [3] one 
concludes that the HSD results for the scaled variances in central A+A colli- 
sions are close to those in N+N collisions. For the SPS energy region all scaled 
variances, u± and u c h, in central A+A collisions are slightly below the N+N 
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Fig. 3. The multiplicities per participant, rii (left), and scaled variances, u>$ (right). 
The solid lines are the HSD results for N+N collisions according to Eq. ([4j) . The full 
circles are the HSD results for central A+A collisions for zero impact parameter, 
6 = 0. The full squares for n_ are the NA49 data [9] for ((ir~) + (K-))/(N P ) in the 
samples of 7% most central Pb+Pb collisions. The full black dots for n c h are the 
PHOBOS data [TO] . The HSD results for uj{ after the subtraction of the contributions 
Acjj ([6]) are shown by open triangles. The dotted lines are the MCE HG model results 
for uji [13] . The HG parameters correspond to the chemical freeze-out conditions 
found from fitting the hadron yields. 

results. The reversed situation is observed for RHIC energies. Thus, the HSD 
results for multiplicity fluctuations are rather similar to those of the WNM. 
This is in agreement with our previous finding [12]. For the samples with a 
fixed number of nucleon participants, Np~ OJ = N l p rg = const, in Pb+Pb colli- 
sions at 158 AGeV, HSD shows fluctuations of the final hadrons close to those 
in N+N collisions at the same energy. This happens to be also valid for most 
central collisions (b = 0) considered in the present study. To estimate the 
influence of participant number fluctuations we calculate the scaled variance, 
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up = Var(Np)/{Np), in the HSD simulations. The participant number fluc- 
tuations are found to be rather small for collisions with 6 = 0. For example, in 
Pb+Pb collisions with b = at 158 AGeV the mean number of participants 
is (Np) = 392, and the scaled variance is up = 0.055 . The additional fluc- 
tuations, AuJi, of ith hadrons due to participant number fluctuations can be 
estimated as [12], 

Aui = rii up . (6) 

The HSD results for u>i after subtraction of the contributions Au^ <^ are 
shown in Fig. [3] by open triangles. The contributions to Ui due to participant 
number fluctuations estimated by Eq. (j^D are small, and they do not explain 
the (positive) difference, Wj(AA) -W;(NN) seen in Fig. [3]at ^/Sjvv = 200 GeV. 

In the statistical model the scaled variances cjj = 1 for the ideal Boltzmann 
gas in the grand canonical ensemble (GCE). The deviations of u>i from unity 
in the hadron-resonance gas (HG) model stem from Bose and Fermi statistics, 
resonance decays, and exactly enforced conservations laws within the canoni- 
cal ensemble (CE) or micro-canonical ensemble (MCE) [T3TH] . Note that the 
statistical model gives no predictions for the energy dependence of hadron mul- 
tiplicities. All yields are proportional to the system volume V which is a free 
model parameter fitted to the multiplicity data at each collision energy. How- 
ever, the statistical model does predict the scaled variances as Ui to become 
independent of the system volume for large systems. In Fig. [3] the scaled vari- 
ances Ui calculated within the MCE HG model along the chemical freeze-out 
line (see Ref. [T3] for details) are presented by the dotted lines: Ui reach their 
asymptotic values at RHIC energies, w±(MCE)= 0.3 and cu c/i (MCE)= 0.6. The 
corresponding results in the GCE and CE are the following: co>±(GCE)= 1.2 
and cj c/i (GCE)= 1.6, cj ± (CE)= 0.8 and lu cJi (CE)= 1.6. The HSD results for u, L 
in central A+A collisions are very different. They remain close to the corre- 
sponding values in p+p collisions and, thus, increase with collision energy as 
Ui oc 7ij. One observes no indication for 'thermalization' of fluctuations in the 
HSD results. This is especially seen for RHIC energies: u;j(HSD)/cG>j(MCE)> 10 
at = 200 GeV. 

Recently first measurements of fluctuations of particle multiplicity [T5fl6lll7] 
and transverse momenta [TB] in A+A collisions have been performed. The 
scaled variance for negatively, positively, and all charged hadrons was mea- 
sured as a function of centrality at SPS [T5fT6"] and RHIC [IT] energies. It 
has been argued in Ref. [12] that the final multiplicity fluctuations seen in 
the NA49 data [15] at the SPS are a straightforward consequence of the fluc- 
tuations in the number of nucleon participants (see also the discussion in 
Ref. [19J). It seems that the same conclusion [20J is valid for the RHIC data 
|17j . In the language of statistical models, the fluctuations of the number of 
nucleon participants correspond to volume fluctuations, hence, they translate 
directly to the final multiplicity fluctuations. To avoid these 'trivial' fluctua- 
tions, one has to select a sample of very central, < 1%, collisions. Such a rigid 
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Fig. 4. Upper panel. The scaled variances uj±° for central Pb+Pb collisions. The 
squares are the NA49 data for 1% most central collisions [21]. The dotted lines 
show the MCE HG model results calculated from full Air scaled variances using 
Eq. ([7]). The full circles present the HSD results in Pb+Pb collisions for 6 = with 
the NA49 experimental acceptance conditions, while the open circles are obtained 
from the 4tt HSD scaled variances using Eq. ([7]). Lower panel. The MCE HG (dotted 
line) and HSD (full circles) results for the 4vt scaled variances u± (the same as in 
Fig. [3]) are shown for SPS energies. 

centrality selection has been recently done for the NA49 data [21] by fixing 
the number of projectile participants, Np~° 3 = A. 

The HG model was compared with the NA49 data [21] for the sample of 1% 
most central collisions at the SPS energies, 20 -r- 158 AGeV in Ref. [13]. It 
was found that the MCE results for u± are very close to the data, they are 
shown by the dashed lines in Fig. HI The NA49 acceptance probabilities for 
positively and negatively charged hadrons are approximately equal, and their 
numerical values are: q = 0.038, 0.063, 0.085, 0.131, 0.163, at the SPS energies 
of 20, 30, 40, 80, 158 AGeV, respectively. In the statistical model the scaled 
variances aj±° for the accepted particles are calculated from u)± in the full 
space according to the acceptance scaling formulae (see Ref. [13] for details): 

= 1 - q + qu>±. (7) 

Note that the energy dependence of lj± c seen in Fig. H] is strongly influenced 
by an increase with energy of the acceptance parameter q: only about 4% of 
the hadrons are detected at 20 AGeV and 16% at 158 AGeV. 
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The comparison of the HSD results for central Pb+Pb collisions (zero impact 
parameter, 6 = 0) with the preliminary NA49 data of 1% most central colli- 
sions, selected by the number of projectile spectators, is presented in Fig. |U It 
demonstrates a good agreement of the HSD results with the preliminary NA49 
data. There are also no essential differences between the MCE HG model and 
the HSD transport model results. Several comments are needed at this point: 
The HSD results for u± c calculated within the NA49 acceptance (full circles 
in Fig. H]) lead to smaller values of uj± c than those obtained with the accep- 
tance scaling formulae (j7j) (open circles in Fig. Hj), i.e. to the violation of the 
acceptance scaling formulae (j7|). This difference may lead to a 10% effect in 
u acc f Qr NA49 acceptance conditions. Thus, the MCE results for uj± c may 
also be about 10% smaller than those obtained from Eq. (|7j) shown in the 
upper panel of Fig. HI 

The lower panel of Fig. [4] demonstrates that the MCE and HSD results for uj± 
at the lowest SPS energy 20 AGeV are accidentally rather close to each other. 
They both are also close to u± in p+p collisions (cf. Figs. [2] and [3]). The HSD 
scaled variances u>i increase with collision energy. In contrast, the MCE Ui 
values remain approximately constant. The ratio of the HSD to MCE values 
of u± reaches about the factor of 2 at the highest SPS energy 158 AGeV. It 
becomes a factor of 10 at the top RHIC energy ^s^m = 200 GeV. However, 
the rigid centrality selection is absent for the available RHIC fluctuation data. 
Due to this reason the participant number fluctuations give a dominant con- 
tribution to Ui. On the other hand, for the SPS data the small values of the 
acceptance, q = 0.04 -=-0.16, and 10% possible ambiguities coming from Eq. fl7|) 
almost mask the difference between the HSD and MCE results (Fig. HI upper 
panel). 

In summary, we have used the HSD transport model to study the particle 
number fluctuations in central A+A collisions at SPS and RHIC energies. 
HSD predicts that the scaled variances Ui in central A+A collisions remain 
close to the corresponding values in p+p collisions and increase with collision 
energy as the multiplicity per participating nucleon, i.e. Ui oc n*. The scaled 
variances u>i calculated within the statistical HG model along the chemical 
freeze-out line show a rather different behavior: uji approach finite values at 
high collision energy. At the top RHIC energy, ^/snn = 200 GeV, the HSD 
values of ^(HSD) is already about 10 times larger than the corresponding 
MCE HG values of c^(MCE) . 

The HSD and HG scaled variances u>i show a different energy dependence such 
that the deviation grows with increasing energy. However, a comparison with 
preliminary NA49 data of very central, < 1%, Pb+Pb collisions at the SPS 
energy range does not distinguish between the HSD and MCE HG results due 
to two reasons: First, the MCE HG and HSD results for cu, at SPS energies are 
not too much different from each other and from uji in p+p collisions. Second, 
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small experimental coverage of the acceptance, q = 0.04 0.16, and 10% 
possible ambiguities coming from the acceptance scaling relation (J7|) make the 
difference between the HSD and MCE HG results almost invisible. 

Thus, our study shows a clear way to distinguish the different models by 
measuring the multiplicity fluctuations. However, the later are very sensitive 
to the experimental acceptance, such that new measurements of Ui for very 
central A+A collisions with large acceptance at both SPS [22] and RHIC 
energies are required to allow for a proper determination of the underlying 
dynamics. 
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